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Abstract. A generalized continuity equation extending the ordinary continuity 
equation has been found using quanternions. It is shown to be compatible with 
Dirac, Schrodinger, Klein-Gordon and diffusion equations. This generalized equation 
is Lorentz invariant. The transport properties of electrons are found to be governed 
by Schrodinger-like equation and not by the diffusion equation. 
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1. Introduction 

The continuity equation governs the conservation of mass/charge/ probability of any 
closed system. This equation involves the spatial distribution of the flux density that 
is related to the temporal variation of the particle density (charge/mass). Ordinarily, 
this equation is derived from the equation of motion. The motion of any continuous 
charge/ mass distribution can be thought of a continuum (field or fluid). The continuity 
equation guarantees that no loss or gain of such quantities. This equation provides us 
with information about the system. The information is carried from one point to another 
by a particle (field) wave. We therefore, trust that such an equation should exhibit this 
property. For electromagnetic field, the information about the motion of the charges 
is carried away by photons. Because electrons have particle-wave nature, the flow of 
electrons is not like the flow of classical objects. The information about the movement 
of electrons is carried by the wave-particle nature that electrons have. The current 
and charge/mass densities are unified if we formulate our basic equations governing 
a particular system in terms of Quaternion. We have recently, shown that Maxwell 
equations can be written as a single quaternionic equation [1] . Moreover, we have found 
that the quaternionic Lorentz force led to the Biot-Savart law. The magnetic field 
produced by the charged particles of the medium is found to be always perpendicular 
to the particle direction of motion. This magnetic field produces a longitudinal wave 
called electroscalar wave propagates with speed of light in vacuum. Moreover, the 
electromagnetic field travels with speed of light in vacuum is the presence of a medium 
if the current and charge densities satisfy the generalized continuity equation. 

We aim in this work at investigating a generalized equation expressing a wave-like 
nature for the charge and current densities, using Quaternions and seek a wave-like 
continuity equation. Such a study will yield the ordinary continuity equation and some 
other additional equations. This system of equations will reveal the wave nature of 
the current and charge (mass/probability) densities. We call the set of these equations 
the generalized continuity equation (GCE). Maxwell, Schrodinger, Klein-Gordon and 
diffusion equations are shown to be compatible with this GCE. Using the GCE, we have 
shown that Ohm's law is equivalent to a Schrodinger-like equation. Hence, the electrical 
properties of metals are consequences of the wave-particle behavior of electrons and not 
merely due to the drift of electrons. 

2. Continuity Equation 

The flow of any continuous medium is governed by the continuity equation. We would 
like here to write the continuity equation in terms of quaternions. The multiplication 
rule for the two quaternion, A = (do , A) and B = (fr , B) is given by [2] 




(1) 
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Accordingly, one can write the quaternion continuity equation in the form 



VJ 




Vpc 2 



V x J 



where 



i d 
cdi 



V 



J = ( ipc , J 



This implies that 



V- J + 



dp 
di 



and 



0, (2) 

(3) 

(4) 
(5) 

(6) 

We call Eqs.(3) - (5) the generalized continuity equation (GCE). Equation (6) states the 
current density J is irrotational. 

In a covariant form, Eqs.(4) - (6) read 

d^ = 0, N yv = d^J v -d v J ll = Q. (7) 

Notice that the tensor N^ v is an antisymmetric tensor. It is evident from Eq.(7) that 
Eqs.(4) - (6) are Lorentz invariant. Now differentiate Eq.(4) partially with respect to 
time and use Eq.(5), we obtain 

Similarly, take the divergence of Eq.(5) and use Eq.(4), we obtain 
1 d 2 J 



- 1 dJ 



V x J = 0. 



(8) 



at 5 



- v v = o, 



(9) 



where p = p(r,t) and J = J(f,t). Therefore, both the current density and charge 
density satisfy the solution of a wave equation propagating with speed of light. We 
remark that the GCE is applicable to any flow whether created by charged particles or 
neutral ones. 



3. Maxwell's Equations 

Maxwell's equations are given by [3] 

o dB 1 OE 

V.£=il, V-B = 0, VxE + — = 0, VxB--— = p J. (10) 

£q Ot C Cft 

Now differentiate Faraday's equation partially with respect to time and employ Ampere's 
equation and the vector identity V x x b\ = V(V • B) — V 2 B to get 

^|_V^„(VXJ). (1!) 
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Similarly, differentiating Ampere's equation partially with respect to time and employing 
Gauss and Faraday equations yields 

1 d 2 E - I (-> 1 dl\ 

7^- VE =-T^" + ^)' < 12) 

According to the GCE the left hand sides of Eqs.(ll) and (12) vanish. This entitles 
us to say that the electromagnetic field travels with speed of light in free space (when 
p = 0, J = 0) and in presence of charge and current; as far as the current and charge 
densities satisfy the GCE. 

Now consider a conducting media defined by the conduction current density (Ohm's 

law) 

J = aE. (13) 
Taking the divergence of Eq.(13) and using the Gauss law, one obtains 

V- J = aV-E = a— . (14) 

Differentiating Eq.(14) partially with respect to time and using Eq.(5) yields 

V p = — — — , — = -D C V p, D c = = . 15 

c 2 e crc at a p a 

Equation (15) is a Schrodinger-like wave equation describing the motion of free electrons. 
A metal can be viewed as a gas of free electrons (Fermions). The thermal properties 
of these electrons are determine by Fermi-Dirac statistics. Thus, Ohm's law satisfies 
Schrodinger equation. This is a very interesting result, since the electron motion is 
governed by Schrodinger equation owing to the wave-particle nature of the electrons. 
Thus, electrons move in a metal as a material wave (de Broglie wave) . 

The Drude has proposed a model of electrical conduction of electrons to explain the 
transport properties of electrons in materials (especially metals) [4]. The model, which 
is an application of kinetic theory, assumes that the microscopic behavior of electrons in 
a solid may be treated classically and looks much like a pinball machine, with a sea of 
constantly jittering electrons bouncing and re-bouncing off heavier, relatively immobile 
positive ions. He has further shown that the Ohm's law is true, where he related the 
current density to the electric field by a linear similar formula. He also connected 
the current as a response to a time-dependent electric field by a complex conductivity. 
This simple classical Drude model provides a very good explanation of DC and AC 
conductivity in metals, the Hall effect, and thermal conductivity (due to electrons) in 
metals. The failure of Drude model to explain specific heat of metal is connected with 
its disregard to the wave nature of electrons. Hence, the motion of electrons in metals 
in not due to diffusion of electrons, the electrons is transported quantum mechanically. 
Hence, a quantum Ohm's law should be considered instead in which the current density 
J is obtained from Schrodinger theory of a free electron (Fermi gas). Thus, quantum 
mechanically, an electron is viewed as a wave traveling through a medium. When the 
wavelength of the electrons is larger than the crystal spacing, the electrons will propagate 
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freely throughout the metal without collision, therefore their scattering result only from 
the imperfections in the crystal lattice of the metal. 

In Schrodinger terminology the diffusion constant D c = i^—-, which shows that the 
diffusion constant is complex. This would imply that the conductivity is complex. In 
Schrodinger paradigm, the current density and charge density are probabilistic quantities 
determined by the wave function of the electron (-0) as [5] 

p = ip*ip, J= -^-(V>*V^-VW) (16) 
2mi 

The current density and the charge here are related by the generalized continuity 
equation. Therefore, the Ohm's current is a probabilistic wave-like current. 

It is thought that electrons move in a metal by drifting, but we have seen here 
the electrical transport properties of metals are propagated by a a material wave rather 
than drifting of electrons. 

Now, apply Eq.(4) in Eq.(13) and employ Gauss law to get 

Let us write the charge density is a separable form 

p(r,t) = p(r)p(t) . (18) 
The time dependence of the charge density p(r, t) is obtained from Eq.(17) as 

pit) = po exp( ) t , po = const. (19) 

The constant r = — is know as the relaxation time, and it is a measure of how fast 
a conducting medium reaches electrostatic equilibrium. The spatial dependence of the 
charge density is obtain from substituting Eq.(17) in Eq.(15). This yields the equation 

V 2 p(r) - -L(r) = 0, X c = £ -^ = ^- = ^. (20) 
A~ a p ca c 

The solution of Eq.(20) is given by 

p(r) = Aexp(-^-) , A = const. (21) 

Therefore, the charge density distribution (space-time) is given by 

p(r, t) = Be'^ (r+Ct) + Fe-^ ( - r+ct) , B,F = const (22) 

This shows that the charges density decays in space and time. However, if uo c is complex, 
we will have an oscillatory charge density solution. This is true if we consider a to be 
complex. The spreading of wavepackets in quantum mechanics is directly related to the 
spreading of probability densities in diffusion (see sec. 4 below). 

Now take the gradient of Eq.(14) and employ Eq.(5) and the vector identity 
V x (V x J) = V(V • J) - V 2 J [6] we obtain 

^ = -D C V 2 J. (23) 
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This equation shows that J is governed by Schrodinger-like equation. Hence, both the 
current density and charged density are governed by a Schrodinger-like equation. 
Now applying Eq.(13) to Eq.(5) and using Eq.(lO), yield 

§ = -^ c V 2 £, (24) 

where we have used the vector identity V x (V x E) = V(V • E) — \7 2 E and the fact 
that V x E = 0. This equation can also be obtained directly from Eq.(23) by dividing 
both sides of Eq.(23) by a and using Eq.(13). 

Equation (20) may define a limiting conductivity of a material (<t ), when we equate 
A c to Compton wavelength of the electron of mass m e . In this case one finds 

(T = — - . (25) 

This amounts to a value of <To = 1.09 x lO 9 ^" 1 ???," 1 . This can be compared with the 
best conductivity of metal (Silver) which is 6.3 x 10 7 f2~ 1 m~ 1 . Thus, the maximum 
possible conductivity is set by quantum mechanics and governed by Eq.(25). Because 
of this the conductivity at zero Kelvin is never infinite but should be limited to a^. 

3.1. Fluid nature of electromagnetic field 

We have recently found that the magnetic field created by the charged particle to be 
given by [1] 

B m = \xE. (26) 
In magnetohydrodynamics this field is govern by the magnetohydrodynamic equation 

— = vx(VxB). (27) 

For any other fluid one has in general 

| = ^{Vx(3), (28) 

where u = V x v is the vorticity of the fluid. This shows that the vorticity of a fluid is 
an integral part to the fluid flow. When an electromagnetic field encounters a charged 
particle, the charged particle induces a vorticity in the field (or fluid). This equally 
occurs in fluids when a fluid encounters an obstacle. I have recently shown that Eq.(6), 
i.e., V x J = 0, implies that u = j% x ( — §|) [7]. Owing to the magnetic field in Eq.(24), 
we have shown recently that this field gives rise to a longitudinal wave (electroscalar 
wave) traveling at speed of light in vacuum besides the electromagnetic fields [1]. 

Hence, it is evident that there is a genuine one-to-one correspondence between 
electrodynamics and hydrodynamics, i.e., they are intimately related. This would 
immediately imply that the electromagnetic fields propagation mimics the fluids flow. 
The electric field resembles the local acceleration and the magnetic fields resembles the 
vorticity of the moving fluid. 
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Using Eq.(13) in Eq.(24), the magnetic field (and vorticity) produced by the 
electrons in a metal vanishes, since B m = = 0, where J = pv. Hence, the Lorentz 
force on conduction electrons is only electric. Hence, no magnetic field inside the 
conductor can be created from the motion of electrons. Therefore, even if we write 
Ohm's law in the general form J = a(E + v x B), we would have obtained the form in 
Eq.(13) too. 



4. Diffusion Equation 

Diffusion is a transport phenomena resulting from random molecular motion of molecules 
from a region of higher concentration to a lower concentration. The result of diffusion 
is a gradual mixing of material. Diffusion is of fundamental importance in physics, 
chemistry, and biology. 

The diffusion equation is given by [8] 

J=-DVp, (29) 

where D is the diffusion constant. This is known as Fick's law. Taking the divergence 
of Eq.(27) and using Eq.(4), one finds 

f t =DV*p. (30) 

This shows that the density p satisfies the Diffusion equation. The normalized solution 
of the Eq.(28) is 

p(xX l = Tmi^-TD? '• (31) 

which when applied in Eq.(27) yields the current density J. It is obvious that 
the current in Eq.(27) satisfies the GCE, viz., Eq.(4) - (6). Differentiation Eq.(27) 
partially with respect to time and employing Eqs.(4) - (5) and the vector identity, 
V(V • J) = V 2 J + V x (V x j) [6] one gets 

^ = ^V 2 J . (32) 
Hence, both p and J satisfy the diffusion equation. 



5. Dirac Equation 

We would like here to show that the generalized continuity equations, viz., Eq.(7) are 
compatible with Dirac equation. To this aim, we apply the current density 4-vector 
according to Dirac formalism, i.e., J M = "07 in Eq.(7). This yields [5] 

d^r - d v ,F = {d^)Y*P + ^Yd^i) - (d'Vh'V - ^(t^V) = 0. (33) 
The first term in the above equation can be written as 

(<9^0) 7 > = (d'V+bVV = (d^+h^V = (7^t/0 + 7 V (34) 
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Hence, 

(<9^) 7 > = @Yd^) + . (35) 

Similarly, 

(d^h'v = (36) 

Substituting these terms in Eq.(31) completes the proof, since for a plane wave 
d^ip = ik^ip. 

6. Klein-Gordon Equation 

It is the equation of motion of scalar particles with integral spin. For such a system the 
current density 4- vector is given by [5] 

2m 

so that the generalized continuity equation becomes 

d,J u - d v J, = ^ [d^d^ - - d^d^ - {d^m = 0, (38) 

and hence, Klein-Gordon equations is compatible with our GCE too. 

7. Schrodinger Equation 

It the equation of motion governing the motion of non-relativistic particles, eg. electrons. 
The current density and probability density in Schrodinger formalism are given by [5] 

p=n>, / = -^*vv> - (w>». (39) 

2mi 

Applying Eq.(6), one gets 
VxJ 



J M =— (37) 



(V x {ip*Vip) - V x (Vip*)ip) . (40) 

2 mi I 

Using the two vector identities V x (fA) = /(V x A) - A x (V/) and V x (V/) = 0, 
Eq.(38) vanishes, and hence one of the GCE is satisfied. We now apply Eq.(5) in Eq.(37) 
to get 

81 -» ft 8 

V(PC2) + dt = V( ^^ + 2m~idi [rV ^ ~ {Vr) ^ ] = °- (41) 
Upon using the Schrodinger equation, 

= ift^t- , i(j*H = -ift^ (42) 
dt dt v 1 

Eq.(39) vanishes for a plane wave equation, i.e., ip(r,t) = Aexp i(k-r — u)t) ,A = const.. 

Therefore, the GCE is satisfied by Dirac, Klein-Gordon and Schrodinger equations. 

This implies that these GCE is indeed fundamental in formulating any field theoretic 

models involving the motion of particles or fluids. This equation will have immense 

consequences when applied to the theory of electrodynamics. Such an application 

will lead to better formulations and understanding of the astrophysical phenomena 

pertaining to the evolution of dense objects. 
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8. Concluding Remarks 

We have derived in this paper the generalized continuity equations and showed how it 
is applicable to the particle flow. We have shown that the the generalized continuity 
equation is lorentz invariant. We have then shown that the basic equations of motion 
are compatible with the generalized continuity equation. The diffusion equation is in 
agreement with the GCE. Moreover, the current and density in Dirac, Klein-Gordon and 
Schrodinger formalisms are compatible with the GCE too. The classical Ohm's law is 
found to be compatible with Schrodinger-like equation. Transport properties of metals 
are shown to be propagated by a material wave rather than drifting and diffusion. 
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